Computational and Structural Biotechnology Journal

Artificial Intelligence (AI) frameworks for automating scientific research have shown strong performance on benchmarks, but their utility for real-world industrial research remains insufficiently characterized. Extending the analysis presented in the first paper of this series, we evaluated the same five advanced AI research frameworks (Kosmos, K-Dense, ToolUniverse, BioAgents from bio.xyz, and the AI Scientist-v2 from Sakana AI) on two more projects of high practical importance for biopharmaceutical development: predicting antibody developability properties with the use of pretrained protein language model embeddings, and modeling non-covalent lipid-RNA interactions in lipid nanoparticles with all-atom molecular dynamics (MD) simulations. The AI frameworks again showed genuine strengths, including unprompted identification of subtle methodological issues, successful use of pretrained protein embeddings, and consistent reporting of p-values and confidence intervals often absent from the original papers. However, no framework approached the scope of the original studies, and severe failures and hallucinations were observed. Our results confirm and extend the conclusion of the first paper that real published research from pharmaceutical companies that we tried to reproduce proved to be considerably harder for current AI frameworks than standard benchmarks suggest.

In vitro fertilization (IVF) laboratories equipped with timelapse incubators generate vast quantities of sequential embryo images, yet the absence of standardized, annotated databases impedes the development of reproducible computational tools for embryo assessment. Here we describe the construction of a standardized time-lapse imaging database comprising 631 normally fertilized zygotes from 218 treatment cycles, integrating timelapse image sequences, patient clinical records, and embryo developmental outcomes. We further present a gradient boosting decision tree (GBDT) ensemble framework that integrates zygote morphokinetic parameters-continuous time-series features extracted via a validated CNN-based segmentation algorithm (US Patent US11210494B2)-with conventional embryo assessment grades (categorical features per the Istanbul consensus). The fusion framework employs equal-weight initialization followed by iterative residual-decreasing training to optimally combine heterogeneous feature types. Ablation analysis demonstrated that the integrated model achieved an AUC of 0.78, significantly outperforming morphokinetics-only (AUC 0.71) and conventional-only (AUC 0.65) models, confirming the complementary value of the two data modalities. The database and fusion framework provide a reproducible foundation for embryo development assessment and are generalizable to other multimodal data integration tasks in reproductive medicine.

Recent advances in artificial intelligence have accelerated the discovery of bioactive peptides by enabling computational exploration of the vast peptide sequence space. However, existing peptide generation approaches generally rely on either distribution-learning models, which generate biologically realistic sequences but do not consistently optimize functional activity, or optimization-based methods, which maximize prediction confidence while often deviating from the underlying distribution of experimentally validated peptides. To address this limitation, a two-phase generative-evolutionary framework is proposed that integrates distribution learning with evolutionary optimization. In the first phase, Variational Autoencoders (VAE), Autoregressive Transformers (ART), and Token Diffusion Transformers (TDT) are used to generate biologically plausible seed peptides. In the second phase, these peptides were used as initial seed for Hill Climbing optimization procedure that iteratively improves fitness function score. The proposed two-phase framework was evaluated using a dataset of experimentally validated IL-2-inducing peptides. Evaluation using independent IL-2 prediction models showed that Autoregressive Transformer combined with Hill Climbing achieved the best overall performance, achieving the mean IL-2 induction confidence score of 0.96 while reducing KL divergence from 2.26 for standalone Hill Climbing to 0.75. A case study on an independent IL-13 inducing peptide dataset showed similar trends, with ART initialized Hill Climbing achieving the mean IL-13 induction score of 0.99 while reducing KL divergence from 1.76 to 0.59. Overall, the framework provides a generalizable approach for balancing functional optimization and distributional realism and can be applied to peptide discovery and data augmentation in imbalanced biological datasets thereby generating high confidence peptides for wet lab validation. HighlightsO_LIProposed a two-phase framework for bioactive peptide generation with potential to address class imbalance in peptide classification tasks. C_LIO_LIPerformed a systematic comparison of distribution-learning and optimization-based approaches for peptide generation. C_LIO_LICombined distribution-learning models for sequence generation with optimization algorithms for improving peptide functional properties. C_LIO_LIDemonstrated the applicability of the proposed framework across multiple bioactive peptide datasets. C_LI

Hydrophobicity is a critical property associated with the risk of non-specific binding, and it is commonly assessed using hydrophobic interaction chromatography retention time. Several computational approaches have been developed to predict antibody developability based on pre-trained language models. Such models can be fine-tuned with limited labeled antibody sequences and, in principle, do not require structural information, which is often challenging to obtain. Nevertheless, few studies have achieved strong performance in hydrophobicity prediction without incorporating structural features. Here, we present a case study of fine-tuning the pre-trained model IgBert to predict antibody hydrophobicity. Using Herceptin as a reference, we performed hydrophobic interaction chromatography retention time experiments and generated Herceptin-adjusted datasets. The fine-tuned model achieved a best R2 of 0.916, underscoring the critical role of rigorous data quality control. We also synthesized and validated 20 commercially available antibody sequences, and the results showed that the predicted hydrophobic properties were correctly reflected. Our findings provide practical guidance and highlight considerations for future applications of fine-tuned pre-trained language models in antibody hydrophobicity prediction. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=189 HEIGHT=200 SRC="FIGDIR/small/742939v1_ufig1.gif" ALT="Figure 1"> View larger version (36K): org.highwire.dtl.DTLVardef@9c814eorg.highwire.dtl.DTLVardef@ed609dorg.highwire.dtl.DTLVardef@62172forg.highwire.dtl.DTLVardef@1e01d37_HPS_FORMAT_FIGEXP M_FIG C_FIG

Oncogenic mutations are known to occur at non-random rates and specific hotspots, yet the structural and thermodynamic factors underlying these hotspots remain poorly understood. This study presents molecular dynamics simulations of 22 cancer hotspots across 10 oncogene families, simulated as histone-DNA complexes. Across 18 of 22 mutations, thermodynamic effects were predominantly localized to within 6 [A] of the DNA-histone interface (mean capture 104%), with van der Waals interactions driving the effect at the thermodynamic extremes, indicating that oncogenic mutations alter nucleosome binding through precise local contact changes rather than global structural rearrangements. Further, analysis revealed a gene-family correlated pattern in thermodynamic stability. RAS family mutations showed a consistent trend toward nucleosome stabilization (mean {Delta}{Delta}G = -4.82 kcal mol-1), while kinase domain mutations trended toward destabilization (mean {Delta}{Delta}G = +40.12 kcal mol-1), a difference reaching statistical significance in this exploratory analysis (Mann-Whitney U, p = 0.008). This interface-specific mechanism, combined with the gene-family-correlated thermodynamic pattern, provides a biophysical framework for understanding nucleosome-level contributions to cancer hotspot mutation biology.

The perception of flavor during food consumption is a complex process governed by the dynamic release of aroma compounds, which is influenced by the interaction between the food matrix and human physiology. While mechanistic models have been developed to simulate these phenomena, their practical application has often been limited by the lack of flexible software tools capable of integrating diverse food structures and real-time physiological data. This study addresses this gap by incorporating in vivo mechanistic models for both liquid foods and melting solid matrices in a modular Python-based framework: digimouth. The modeling approach treats the oral, pharyngeal, and nasal cavities as a series of interconnected reactors with mass transfer governed by ordinary differential equations (ODEs). These models explicitly integrate key physiological parameters, including salivary flow, bidirectional respiratory airflow, and discrete events such as mastication and swallowing. For solid foods, the framework incorporates additional mechanisms for matrix dissolution, temperature-dependent melting, and dynamic surface area change. The models were validated using experimental data from 30 subjects who consumed aroma solutions and flavored gelatin gummy discs while monitored by PTR-MS (Proton Transfer Reaction - Mass Spectrometry). The simulations successfully reproduced the main kinetic features and release profiles observed experimentally for both types of matrices. Despite some limitations related to the self-reporting of physiological events and the non-inclusion of initial orthonasal mass transfer in the models, the developed software package proved to be a robust and evolutive tool. It offers a foundation for exploring inter-individual variability and the impact of complex food structures on sensory perception.

Artificial intelligence (AI) is increasingly permeating the drug development pipeline. Numerous algorithms for accelerating this multi-stage and multi-task process have been constructed, which depends heavily on expert design and labor-intensive task-specific optimization. Given that AI-driven acceleration of drug development is recognized as a cumulative, often synergistic, effect across multiple stages, the autonomous evolution of existing algorithms across the entire pipeline is demanded to achieve a holistic advancement. Here, we present DrugEvolve, a multi-role large language model system for systematic and autonomous algorithm evolution in drug development. DrugEvolve realizes a closed-loop evolution process by incorporating Researcher, Engineer, and Analyst domains, and enables an iterative design, implementation, evaluation, and refinement of algorithm by leveraging scientific knowledge and accumulated evolutionary experience. Across eleven representative tasks spanning target identification, drug discovery, preclinical study, and clinical trial, DrugEvolve autonomously evolved the corresponding task-specific algorithms and achieved substantial performance enhancement on 120 benchmark test sets. Moreover, it showed robust generalizabilities across heterogeneous data modalities (ranging from biological sequence and graph to molecular topology and textual language), and realized gains in both predictive and generative tasks. Collectively, this AI system can serve not only as an algorithmic infrastructure for drug development, but also as a transferable paradigm for broader scientific domains.

Alzheimer's disease is characterized by the accumulation and aggregation of amyloid-{beta}(A{beta}), but the molecular mechanisms linking environmental and infectious factors to A$\beta$ conformational changes remain incompletely understood. Herpes simplex virus type 1 (HSV-1) has been proposed as a potential contributor to AD pathology, and interactions between the viral glycoprotein B (gB) and A$\beta$ may influence the conformational behaviour of the peptide. Molecular dynamics (MD) simulations provide atomic-scale information on such interactions, but conventional structural descriptors may not fully capture changes in the organization of residue interaction networks. Here, we introduce a graph-geometric framework based on Forman-Ricci curvature to characterize the evolution of residue interaction networks during MD simulations. Each simulation frame is represented as a residue interaction graph based on C--C contacts, and residue-wise curvature profiles are analysed across time. We apply the framework to A{beta}1-42 in isolation and in complex with HSV-1 gB. Conventional MD analyses indicate stable association of the simulated complex, favourable interaction energetics, and conformational changes in A{beta}, including a transition from -helical structure toward {beta}-turn-rich conformations over the simulated timescale. Forman-Ricci curvature reveals pronounced and spatially localized remodelling of the A{beta} residue interaction network in the complex, with the strongest changes concentrated in the C-terminal region. These regions also exhibit reduced temporal curvature fluctuations and progressively distinct geometric behaviour throughout the simulation. Hierarchical clustering further identifies cooperative groups of residues with coordinated curvature dynamics, including a prominent C-terminal domain. Together, these results demonstrate that Forman-Ricci curvature provides a complementary description of biomolecular dynamics by capturing changes in the geometric organization of residue interaction networks that are not directly represented by conventional structural descriptors. The framework provides a general computational approach for studying network-level structural remodelling in protein molecular dynamics and offers a quantitative perspective on the conformational consequences of HSV-1 gB--A{beta} association.

Antibodies generally exhibit high specificity for their cognate epitopes, but structural and physicochemical similarities between distinct epitopes can enable an antibody to recognize different antigens, resulting in cross-reactivity. This property can be exploited for antibody repurposing. To identify epitopes that share such similarities, both sequence- and structure-based approaches can be employed. In this context, 3D Zernike descriptors provide a compact representation of protein surface geometry as numerical feature vectors, enabling quantitative comparisons independently of structural alignment and orientation. Thus, this study aimed to evaluate the application of 3D Zernike descriptors for the structural clustering of antibodies and epitopes and to explore their use in antibody repurposing for the recognition of new targets. To this end, antibody binding sites previously associated with recognition of similar epitopes were analyzed at different structural levels, considering the CDRs, CDRH3, and complete paratopes. Surface similarity was subsequently quantified by calculating the Euclidean distance between their corresponding 3D Zernike feature vectors. Performance was benchmarked against SPACE2. Additionally, different distance thresholds were evaluated based on their ability to recover antibody pairs recognizing the same epitope. The paratope-based approach provided the best balance between the number of identified pairs and precision at a distance threshold of 2.7, whereas epitope clustering showed robust performance up to a distance of 3.0. At these thresholds, the 3D Zernike descriptors identified a greater number of functional pairs than SPACE2 while maintaining comparable precision and identifying complementary sets of antibody pairs.. BTaken together, these findings support the use of 3D Zernike descriptors for structural clustering of antibodies and epitopes and for guiding antibody repurposing G, a highly lethal zoonotic pathogen. Structural screening identified three antibodies with epitopes similar to the NiV target that also showed a consistent binding preference for the target epitope in molecular docking assays. Notably, one candidate, originally directed against a SARS-CoV-2 epitope, formed a stable complex with the NiV epitope, remaining within the 5 [A] RMSD threshold during heated molecular dynamics simulations and emerging as a potential cross-reactive candidate.These results support the use of this computational framework for biopharmaceutical discovery against emerging targets. Taken together, these findings support the use of 3D Zernike descriptors for structural clustering of antibodies and epitopes and for guiding antibody repurposing.

The attention mechanism offers the possibility for data-driven discovery of biological principles. However, for important protein families such as human olfactory receptors, the extent to which attention can associate with biologically meaningful key regions lacks systematic validation. In this study, using human olfactory receptors (ORs) as a model, we constructed CrossVOI, a VOC-OR interaction prediction framework based on protein language models and cross-attention, achieving predictive performance superior to existing methods. Furthermore, we systematically analyzed the attention distributions of CrossVOI and found that attention not only focused on ligand-binding interfaces and evolutionarily conserved sites, but also to some extent identified certain dynamically regulated regions. In summary, we propose CrossVOI, currently the best-performing framework for VOC-OR interaction prediction, and analyze the interpretability of the attention mechanism for human ORs. This study provides insights into the interpretability of protein function prediction methods and is expected to contribute to the exploration of attention mechanisms in biological mechanisms, and provide assistance for large-scale screening and mechanistic analysis of olfactory receptors.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterised by progressive motor neuron degeneration. Mutations in the SOD1 gene represent the second most common genetic cause of ALS (ALS), and distinct SOD1 missense variants present with markedly different clinical profiles. A4V leads to an aggressive form of the disease (median survival [~]1y), H46R confers a mild, slowly progressive course and I113T exhibits an intermediate phenotype. The molecular basis by which these mutations produce divergent clinical outcomes remains poorly understood. We performed extensive classical molecular dynamics simulations of wild-type SOD1 and the three ALS-associated variants in the apo monomeric state to attempt to investigate the mechanisms behind such phenotypic differences. Structural stability, global compactness, and conformational flexibility, as well as analysis of collective motions between residues and estimation of free energy, were assessed. The H46R, A4V, and I113T variants exhibited distinct dynamic behaviours, highlighting differences in structural stability, local flexibility, and intramolecular interactions. These findings suggest that specific structural regions may contribute differently to protein dysfunction and could represent key elements for understanding the relationship between molecular dynamic properties and the differing clinical severity associated with these variants. Most strikingly, H46R exhibited exceptional structural stability across every analytical level, the lowest global deviation, most attenuated local flexibility, strongest internal dynamic coordination, and the deepest, most confined free energy basins of any system examined. This convergent multi-layered evidence of structural restraint provides a compelling mechanistic basis for the mild and slowly progressive clinical course of H46R ALS, suggesting that enhanced conformational rigidity, rather than bulk destabilisation, is the defining biophysical feature of this variant, and that its pathogenic mechanism operates through a route fundamentally decoupled from the aggregation-driven toxicity that characterises the more aggressive SOD1-ALS mutations.

Antibodies provide programmable molecular recognition, whereas enzymes enable repeated chemical transformation. Catalytic antibodies seek to combine these properties within a single protein scaffold. However, conventional approaches based on transition state analogue immunisation, library screening or local mutagenesis provide limited control over the atomic arrangement of catalytic residues. They also frequently produce antibodies that bind substrates without supporting efficient chemical turnover. Recent advances in generative protein design have enabled the construction of antibody complementarity determining regions and the scaffolding of functional motifs under structural constraints. A systematic strategy for transferring experimentally supported enzyme active site geometry into antibody variable domains is still lacking. Here, we present a computational framework that treats antibody and enzyme structures as distinct but complementary inputs. Developable Fv or VHH structures provide the immunoglobulin scaffold. Enzyme complexes containing substrates, products or transition state analogues provide catalytic residues, ligand conformations, metals, cofactors and key water networks. The selected catalytic atoms are mapped into antibody complementarity determining regions, while the surrounding loops are reconstructed using antibody compatible representations and constrained all atom diffusion. Sequence design and structural back prediction are followed by filters for antibody folding, catalytic geometry, ligand positioning, conformational stability and developability. The framework avoids direct fusion of intact enzymes and antibodies. Instead, it transfers only the local geometry required for catalysis. This separation of scaffold selection from catalytic motif selection creates a testable route for determining whether natural enzyme chemistry can be embedded within antibody formats. It also provides a practical basis for evaluating substrate binding, chemical conversion, product release and catalytic turnover as separate design objectives. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=91 SRC="FIGDIR/small/740676v1_ufig1.gif" ALT="Figure 1"> View larger version (45K): org.highwire.dtl.DTLVardef@16ef6d6org.highwire.dtl.DTLVardef@f7c23org.highwire.dtl.DTLVardef@9ee50borg.highwire.dtl.DTLVardef@1cf42d4_HPS_FORMAT_FIGEXP M_FIG C_FIG

Background: Peanuts are among the most prevalent food allergens responsible for anaphylaxis, particularly in children. Allergic responses may be mitigated by disrupting molecular interactions between immunoglobulin-E (IgE) and its binding partners: peanut allergen (Ara h 2) and IgE receptors (Fc{epsilon}RI). However, the structural dynamics of IgE upon engaging these factors are not yet fully elucidated. Objective: To characterize how Ara h 2 and/or Fc{epsilon}RI binding influence IgE flexibility and its interdomain communications. Methods: The structural dynamics of full-length IgE in various unbound and bound states were examined by combining molecular dynamics (MD) simulations with cross linking mass spectrometry (XLMS). Causal relationship between IgE domains were characterized to infer allosteric communication pathways, which were assessed through hydrogen-deuterium exchange mass spectrometry (HDX-MS). Results: Non canonical bent IgE conformations were identified. We demonstrate that Fc{epsilon}RI binding immobilizes IgE predominantly by stabilizing and restricting Fc flexibility, whereas Ara h 2 binding induces specific conformational changes within the Fc{epsilon}RI-binding region. Our results further suggest that the C{epsilon} domains retain sufficient freedom to enable causal dynamics by the Fabs in a binding-dependent manner, thereby modulating the receptor engagement potential of IgE. Notably, concurrent binding of Ara h 2 and Fc{varepsilon}RI to IgE shifts the allosteric communication between the Fabs and Fc regions from the C{epsilon}3 to C{epsilon}4 domain. Conclusion: Our findings provide new insights into the dynamic regulation of IgE and its role in allergic responses. This could serve as a structural framework for understanding IgE signaling in peanut allergy and suggest potential targets for therapeutic interventions

AlphaFold3 has shown promise as a tool for predicting antibody-antigen binding, yet its performance across large datasets has not been fully characterized. In this study, 3401 experimentally validated antibody-antigen complexes were sourced from the Structural Antibody Database and screened alongside 23798 negative controls to benchmark AlphaFold3s binding prediction capabilities. Confidence metrics including Predicted Aligned Error and Interface Predicted Template Modeling score were used to achieving a maximum recall of 53% at 100 inference seeds. Several factors were found to influence prediction accuracy: a notable bias was observed toward antibodies derived from X-ray crystallography structures versus those from electron microscopy, and positive prediction rates were found to decrease with increasing target protein size and surface area. In contrast, neither the amino acid composition or lengths of the complementarity determining regions, nor training data leakage were found to introduce significant bias. An innate false positive rate of approximately 3% was identified, with AF3 shown to hallucinate plausible binding interfaces across the surface of decoy targets while avoiding disordered regions. Epitope mapping using DockQ, epitope shift, and antibody displacement revealed that approximately 34% of false negatives retained the correct epitope location despite poor structural alignment, suggesting that conformation refinement tools could recover additional true binding predictions. These findings provide a comprehensive characterization of AlphaFold3s strengths and limitations for antibody screening in computational drug discovery. Key MessagesO_LIAlphaFold3 has a recall of 50% and an innate false positive prediction rate of 3%. C_LIO_LIFalse negative predictions can still feature the correct epitope despite poor RMSD. C_LIO_LIFactors such as disorder and target size impact accuracy. C_LI

The modeling of aroma release in food science is a complex, multidisciplinary field integrating mass transfer theory and principles of physical chemistry. While mechanistic models have long been used to understand volatile release, their practical adoption has been hindered because existing implementations are often private, "in-house" tools that lack reproducibility and extensibility. This paper introduces digimouth, an open-source Python framework designed to provide a unified, modular, and reproducible environment for simulating and fitting aroma release kinetics. A key innovation of digimouth is its handling of events and stages. Events represent discrete, instantaneous changes in the system (such as opening a bottle or swallowing), while stages allow for multi-phase simulations where experimental conditions or physical dynamics change over time. To demonstrate the tool's utility, the authors implemented established in vitro convection models and performed sensitivity analyses on parameters such as air-water partition coefficients and mass transfer coefficients. The package also includes an optimization module based on the least-squares method to fit parameters to experimental data. Validation against real-time data from Proton Transfer Reaction-Time of Flight- Mass Spectrometry (PTR-ToF-MS) shows a satisfactory fit, with an average error of approximately 2%. Ultimately, digimouth facilitates robust model development and hypothesis testing, offering the food science community a collaborative tool to better understand the mechanisms of flavor release.

Codon optimization uses synonymous sequence changes to improve the expression and therapeutic performance of nucleic acid-based medicines. Masked language models (MLMs) have recently been proposed as alternatives to traditional, frequency-based codon optimization approaches, yet whether they offer a meaningful advantage over such simpler methods remains unclear. Here we benchmark three prominent MLMs - CaLM, EnCodon and CodonTransformer - across backtranslation fidelity, sequence generation and nine molecular phenotype prediction tasks, and experimentally evaluate model-designed sequences using a secreted embryonic alkaline phosphatase (SEAP) reporter. The models differed markedly in amino-acid fidelity and generated distinct synonymous sequence variants. However, no single model performed best across all benchmark tasks and simple sequence features remained competitive in several settings. Our interpretability analysis revealed that the models integrate a large window of codon context for making predictions, as opposed to frequency-based approaches. Our in vitro data showed that MLM-designed variants outperformed conventional and commercial-vendor-derived sequences in both transient and stably integrated expression, supporting the models ability to capture translational context beyond codon frequency. Together, our results establish MLMs as effective and complementary tools for codon optimization and suggest that sampling across multiple models may improve the likelihood of identifying high-performing therapeutic sequences.

T cell receptor (TCR) specificity prediction is critical for understanding adaptive immunity and for enabling the development of therapeutic interventions. Current machine learning models for predicting TCR-peptide-MHC interactions perform poorly, despite large efforts by many groups and advances in protein language models. This work identifies and addresses the fundamental issues that might be hindering progress in this field, focusing on current approaches for constructing training and validation sets, particularly on the strategies used by commonly used methods, which generate negative data either by shuffling TCRs between epitopes or by pairing epitopes with TCRs from background datasets of healthy donors. Here, we visualize and quantify the overlap between positive and generated negative TCR sets produced by both strategies using a range of sequence embedding approaches, including simple encodings (one-hot and BLOSUM62), general protein language models (ESM-2 and ProtGPT2), and a TCR-specific embedding model (TCR-BERT). Using paired {beta}TCR data from VDJdb for two well-studied epitopes (GILGFVFTL and KLGGALQAK), we train Random Forest classifiers on each embedding representation and evaluate their ability to distinguish positive from negative data examples. We find that general sequence embedding methods struggle to differentiate between positive and negative TCRs. We further show that the choice of negative data generation method significantly impacts classifier performance. Our findings highlight fundamental limitations in current TCR specificity training data and general protein language models, and underscore the need for TCR-specific embeddings, experimentally validated negative datasets, transparent reporting of negative data construction, and epitope-specific modeling approaches.

Abstract Background: Immune-mediated skin diseases frequently co-occur, but the molecular architecture connecting clinically distinct disorders, and the extent to which molecular proximity predicts shared therapeutic opportunity, remain incompletely defined. We constructed a multilayer skin diseasome integrating transcriptomic overlap, pathway enrichment, protein-interactome topology, curated disease associations, real-world clinical co-occurrence, and drug-perturbation signatures across psoriasis (PS), atopic dermatitis (AD), hidradenitis suppurativa (HS), vitiligo (VL), and chronic spontaneous urticaria (CSU), using sodium lauryl sulfate-induced irritant contact dermatitis (SLS-ICD) as a mechanistically distinct barrier-injury comparator. Methods: Differentially expressed genes (DEGs) were derived from publicly available NCBI Gene Expression Omnibus datasets and compared using exclusive pairwise and higher-order intersections. Gene Ontology biological-process enrichment characterized shared and condition-specific programs. Interactome-based network separation (S_AB) quantified the relative topological proximity of disease modules. DisGeNET enrichment generated candidate disease +associations, a selected subset of which was examined descriptively for patient co-occurrence in TriNetX. Candidate compounds were ranked by reversal of disease-associated expression signatures against LINCS L1000 perturbation profiles. Results: The largest exclusive DEG intersection occurred between PS and HS (2,060 genes), followed by AD-PS-HS (490), VL-PS (464), AD-PS (402), and VL-PS-HS (379). Fourteen genes were shared across PS, AD, HS, VL, and CSU but not SLS-ICD. Neutrophil-associated processes recurred across the immune-mediated conditions, whereas SLS-ICD was dominated by keratinization and epidermal-differentiation programs. PS and HS had the lowest observed network-separation value (S_AB = -0.594); SLS-ICD had positive S_AB values in all five comparisons, although this separation did not persist once module sizes were harmonized. DisGeNET linked the index conditions to diverse immune, cardiovascular, neurologic, and neoplastic categories; 42 of these candidate diseases were observed to co-occur with an index condition in the queried TriNetX cohorts. Drug-signature overlap was greatest between PS and CSU, despite their marginal interactome separation, indicating that interactome proximity and predicted pharmacologic reversal capture different dimensions of disease relatedness. mTOR/PI3K-pathway inhibitors were prominent among the psoriasis reversal candidates, while methotrexate and mycophenolate mofetil emerged from both the CSU signature and the curated gene set associated with HIV infection. Conclusion: Multilayer integration identified a connected molecular core among five immune-mediated skin diseases and a distinct SLS-ICD injury module. The discordance between interactome proximity and drug-signature overlap suggests that no single layer adequately defines disease relatedness. These findings provide a hypothesis-generating framework for prioritizing disease links and repurposing candidates, but controlled epidemiologic analyses and experimental validation are required before clinical or therapeutic inference.

Changes in pH play an important role in human physiology and cell behaviours, evidenced by studies of tissues and gene expression, and the biophysics of individual systems. Bridging between these scales to establish pH-sensing networks has only been partially explored. Here, amino acid conservation and structural context, and benchmarking with known pH-sensors, were used to predict histidine residues within the human proteome that could modulate pH-dependence. Proteins containing these histidines are enriched in many categories, most notably nucleic acid-binding and transcriptional regulation. DNA phosphate backbone interactions are proposed to play a major role in modulating some pH-sensors. Transcription factors with potential pH-sensors and the largest number of documented target genes are highlighted, including those involved in circadian and developmental processes. Certain CDK-cyclin pairs also feature, along with components of SWI/SNF and cohesin complexes. Focussing on histidines, including those that are solvent accessible in AlphaFold protomer models, this study generates benchmarked predictions for possible pH-sensors in the human proteome.

Motivation: MicroRNA-target interaction prediction remains challenging because many existing tools provide prediction scores or ranked candidate lists without making the supporting evidence easy to interpret or relate to a specific biological context. Results: Here, we developed miRAssist, a context-aware evidence-integration framework for interpretable miRNA-target prioritization. miRAssist integrates six evidence families, including sequence complementarity, thermodynamic stability, sequence conservation, target-site accessibility, functional binding, and functional repression. A sequence-defined candidate universe was generated, resulting in 280,917 candidate interactions. Using miRTarBase-supported interactions as known-positive labels, six supervised scoring approaches were evaluated using a grouped train/test split by miRNA. Random forest showed the strongest performance and was selected. miRAssist also produced stronger known-positive enrichment than established miRNA-target prediction models in the evaluated benchmark. An LLM-assisted interface further supports natural-language database querying and evidence-grounded summarization of prioritized candidates.